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abstract 


Since  1946,  the  Defense  Mapping  Agency  (DMA)  has  been  participating 
with  Mexico  and  each  of  the  Central  American  countries,  as  a  result 
of  the  collaborative  program  under  the  Inter  American  Geodetic 
Survey  for  the  production  of  cartographic,  geodetic,  and  geophys¬ 
ical  products*  The  geodetic  leveling  network,  which  developed  from 
these  mapping  agreements  with  each  country,  now  contains  over 
50,000  kilometers  of  leveling.  Gravity  observations  have  also 
been  observed  over  a  large  portion  of  these  level  lines.  Most  of 
these  leveling  data  were  adjusted  by  DMA  in  regional  blocks  as 
the  basic  leveling  networks  were  completed.  Each  regional  block 
was  adjusted,  holding  mean-sea-level  at  the  tidal  stations  equal 
to  zero.  Elevations  of  previously  adjusted  blocks  were  held 
fixed  in  each  successive  block  adjustment.  A  test  adjustment  (i960) 
combining  all  of  the  first-order  level  links  from  the  1959  to 
1972  regional  adjustments  was  performed  for  analysis.  Results  of 
this  1980  Special  Adjustment  along  with  the  current  status  of 
the  leveling  data  in  DMA  files  are  presented. 


1 


INTRODUCTION 


Survey  records  for  some  50,000  kilometers  of  geodetic  leveling 
in  Mexico  and  Central  America  are  available  at  the  Defense  Mapping 
Agency  Hydrographic/Topographic  Center  (DMAHTC).  These  data  were 
observed  in  the  1946-72  period  and  are  the  result  of  cooperative 
mapping  agreements  with  each  country.  Most  of  these  leveling  data 
were  observed  using  first-order  specifications  by  both  personnel 
from  the  Inter  American  Geodetic  Survey  (IAGS)  and  from  the  country 
involved  (  Skaggs  197 9  ). 

2.  THE  1959  THROUGH  1972  ADJUSTMENTS 

As  the  field  observations  were  completed,  they  v/ere  furnished 
to  DMA  through  IAGS,  where  the  data  were  preprocessed  and  checked 
against  accuracy  specifications.  Double  run  leveling  was  used  for 
first-order  work.  The  allowable  discrepancy  between  forward  and  -i 

backward  leveling  between  bench  marks  is  4mm •/l’,  where  L  is  the 
distance  in  kilometers  between  bench  marks.  At  DMA,  the  latitude  for  J 
each  bench  mark  was  scaled  from  available  map  sources,  and  the 
observed  leveling  observations  corrected  for  the  effect  of  non- 
parallelism  of  equipotential  surfaces  (orthometric  correction) 
based  on  normal  gravity. 

* 

The  leveling  data  were  combined  in  regional  blocks  and  adjusted 
by  least  squares  as  the  data  were  received.  These  adjusted  regional 
networks  are  shown  in  figure  1.  The  MX  1959  block  (covering  Mexico) 
and  the  PN  1959  block  (covering  Panama)  were  adjusted  in  the  year 
1959.  Each  new  block  was  adjusted  to  fit  the  previously  adjusted 
blocks.  Most  of  the  tidal  stations  along  both  coastlines  were  held 
fixed,  continuing  the  method  used  in  the  adjustment  of  the  leveling 
data  in  the  United  States  and  Canada  in  1929*  Although  these  eleva¬ 
tions  are  adequate  for  controlling  trigonometric  leveling  networks 
and  large  scale  mapping  projects  they  are  not  accurate  enough  for 
all  modern  engineering  and  scientific  studies.  Some  level  lines  are 
warped  by  holding  mean-sea-level  at  the  tide  gauges  equal  to  zero. 

Thus  the  accuracy  of  leveling,  observed  using  first-order  procedures 
is  downgraded  by  the  adjustment  procedure.  An  example  of  the  networks 


distortion  can  be  seen,  when  we  compare  free  adjustments  verses 
adjustments  where  the  mean-sea-level  was  held  fixed  at  zero  at 
the  tidal  stations. 


REGIONAL 

BLOCK 


CA  I960 
CA  1967 


KILOMETERS 

6697 

4766 


STANDARD  ERROR  OF  UNIT  i.  EIGHT 
FREE  ADJUSTMENT  MSL  FIXED 

2.6  4.3 

2.5  3.5 


The  standard  error  of  unit  weight  is  calculated  by  £q(3). 


3.  LEVELING  DATA  RECEIVED  AFTER  1972 

By  .1972  most  of  the  basic  leveling  networks  had  been  completed. 
Regional  block  adjustments  of  these  data  by  DMA  had  also  been  com¬ 
pleted.  Survey  efforts  by  each  country  since  1972  has  been  in  the 
form  of  extension  and  densification  surveys.  Re-leveling  surveys 
have  also  been  performed  to  keep  the  networks  current.  DMA  continues 
to  process  leveling  data  only  on  a  case-by-case  basis  since  most 
of  the  countries  have  become  self-sufficient.  Approximately  10,000 
kilometers  of  leveling  data  have  been  received  at  DMA  since  1972 
and  these  data  have  not  been  adjusted.  These  level  lines  are 
shown  as  broken  lines  in  figure  1.  A  review  of  progress  charts 
indicate  another  7,000  kilometers  of  leveling  data  as  being  com¬ 
plete  or  near  completion.  There  lines  are  shown  as  dotted  lines 
in  figure  1. 

DMA*s  largest  current  effort  in  level  processing  is  along 
the  Vvest  Coast  of  South  America,  where  the  previous  processed 
regional  blocks  are  being  combined  into  a  single  network  and 
adjusted.  New  survey  data  are  being  processed  and  incorporated 
into  this  combined  network.  Preliminary  results  of  this  work  have 
been  published  (  Etc ty  1977  ). 

4.  THE  1980  SPECIAL  ADJUSTMENT 

To  provide  a  current  evaluation  of  the  leveling  in  Mexico  and 
Central  America,  the  previously  processed  level  links  from  the 
regional  adjustments  were  keypunched  and  assembled  into  a  single 
simultaneous  network  adjustment.  The  network  is  continuous  from 
5  ties  at  the  U.S. -Mexico  border  through  Central  America  to  Panama. 


The  network  was  adjusted  by  DMA  using  the  computer  program 
LEVEL  (  Googz  S  LeAoy  1970  ) .  Weights  were  assigned  to  the  link 

observations  based  on  distance.  The  a* priori  standard  error  for 
each  link  is  played  by  the  square  root  of  the  distance; 

0-i  -Vl^  ,  Lt  i  s  the  link  distance  in  km.  (1) 

and  the  observation  weight: 

V  1  /cl  -  1  /  t,  (2) 

The  normal  equations  were  inverted  to  provide  error  analysis 
upon  completion  of  the  least  squares  solution.  The  standard  error 
of  unit  weight  is  calculated  by  the  formula:: 

1V?V 


<*o' 


(3) 


m-n 


where: 

m  »  number  of  observations  (links) 
n  «  number  of  unknowns  (elevations) 

P  »  weight  of  observation 
V  -  least  squares  residual 

First,  the  level  links  were  combined  into  circuits  and  each  loop 
closure  checked  against  the  4mmVL  criterion,  Six  of  the  72 
circuits  had  closures  over  the  allowable,  however  only  one  was 
significant.  The  circuit  length  is  274  kilometers  with  an  allowable 
error  of  66mm  and  a  closure  of  136mm.  The  reason  for  this  descrep- 
ancy  has  not  been  found.  A  very  large  accumulation  of  systematic 
error  is  a  possibility.  The  circuit  also  falls  in  an  area  of  high 
seismic  activity  and  land  mass  shifts  during  the  leveling  period  is 
also  possible.  The  line  (E64  to  E93),  shown  in  figure  2  by  the 
symbol  was  not  used  in  the  least  squares  adjustment.  The  average 
circuit  distance  in  the  network  is  326  km.  The  average  misclosure 
of  circuits  is  1.9mm-</P,  where  P  is  the  perimeter  distance  for  the 
circuit. 

A  least  squares  adjustment  of  the  links  was  performed  holding 
only  one  station  fixed,  Bench  Mark  A  630  at  the  U.S.  -  Mexico 
border.  The  standard  error  of  unit  weight  from  the  adjustment  was 
2.3rcnVT.  A  statistical  analysis  of  the  residuals  from  the  adjust- 


ment  are  given  in  Table  1.  Nine  normalized  residuals  out  of  the 
246  links  were  flagged  as  falling  over  3.0.  The  normalized  re¬ 
siduals  (vVP)  are  plotted  against  the  Gaussian  normal  distri¬ 
bution  curve  for  convenience.  The  actual  distribution  curve  is 
shown  in  Table  1 .  It  departs  from  the  normal  curve  in  the  0.2  to 
0.3  increment  range.  The  average  link  correction  from  the  least 
squares  adjustment  was  9mm.  As  seen  in  figure  2,  there  is  a 
large  difference  in  the  length  of  each  link.  The  average  link 
correction  tvhen  related  to  the  distance  was  1  .OVL*,  L  is  the 
link  distance  in  kilometers. 

The  standard  error  of  the  computed  elevations  from  the  free 
adjustment  propagate  from  zero  at  Bench  Mark  A  630  in  the  United 
States,  which  ’was  held  fixed,  to  a  maximum  of  121mm  at  bench  marks  ■'* 
Eastern  Panama.  These  adjusted  differences  in  elevation  were  used 
together  with  the  spur  line  differences  to  the  tidal  stations,  to 
determine  the  relative  heights  of  local  mean- sea-level  along  both 
coastlines.  After  the  adjustment  the  elevations  were  all  uniformly  j 
shifted  to  Port  Isabel,  Texas  so  it  would  be  -0.23  meters  and  j 

correspond  to  the  1963  Special  Adjustment  by  NGS  (Bftaaizn  S  McCamks  79fi 
The  local  me an- sea- level  variations,  indicated  by  this  1930  Special  j 
Adjustment,  and  selected  relative  errors  between  the  tidal  bench  ] 
marks  are  given  in  figure  2.  Notice  that  the  mean- sea- level  on  the  j 
Pacific  seems  to  be  about  2  decimeters  higher  than  mean-sea-level  j 
on  the  Gulf  and  Caribbean  Sea  side.  In  the  United  States  the  mean-  j 
sea-level  on  the  Pacific  Coast  was  found  to  be  about  6  decimeters  j 
higher  than  mean-sea-level  on  the  Atlantic  Coast.  j 

Adjustments  of  the  overall  net,  by  constraining  the  tidal  : 

bench  marks  ’with  varing  weights  (based  on  the  tidal  observation  i 

periods),  resulted  in  adjustments  with  the  standard  error  of  unit  ] 

weight  ranging  from  3*9  to  5*6.  This  increase  in  Oi  is  caused  by 
the  difference  in  height  of  mean- sea-level  at  the  two  Coasts.  The 
network  is  forced  to  tilt  by  these  constraints.  A  better  solution 
would  be  to  constrain  the  differences  in  elevations  between  tidal 
bench  marks.  j 


.  THE  NORTH  AMERICAN  REDEFINITION 

The  need  for  a  redefinition  of  the  vertical  reference 
system  in  Mexico  and  Central  America  is  desirable.  The  present 
vertical  system,  referenced  to  mean- sea-level  will  not  meet  ail  0f 
the  requirements  of  future  surveying  and  scientific  needs.  New  lev 
ing  observations  will  have  difficulty  in  fitting  the  previously 
determined  elevations.  Almost  as  important  as  the  method  of 
adjustment  i3  the  upgrading  of  the  network,  which  can  be  ac¬ 
complished  by  applying  refined  leveling  corrections  that  have 
not  previously  been  applied  {  Holdahl  1979  )•  The  refraction 
correction  should  be  applied  to  correct  for  the  bending  of  the 
optical  ray  path.  The  astronomic  correction  can  be  applied  to 
account  for  the  tilting  of  the  equipotential  surfaces,  relative 
to  the  terrain,  caused  by  the  changing  positions  of  the  sun  and 
moon.  The  effect  of  real  gravity  has  also  not  been  taken  into 
account  in  the  previous  adjustments.  To  correct  for  the  effect  of 
non-parallelism  of  equipotential  surfaces,  the  crthometric  cor¬ 
rection  was  applied  based  on  normal  gravity.  Observed  gravity 
data  does  exist  for  most  of  the  leveling  in  Mexico  and  (Central 
America.  The  significance  of  these  refinement  corrections,  as 
they  would  apply  to  Mexico  and  Central  America  network,  is  under 
study  at  DMA. 

6.  CONCLUSIONS 

The  basic  leveling  network  in  Mexico  and  Central  America 
was  observed  using  first-order  field  procedures.  The  network  is 
continuous  from  the  Mexico- United  States  border  to  Panama. 

Although  the  results  of  this  19<50  Special  Adjustment  show  some 
departure  from  expected  first-order  results,  the  net  could  be 
upgraded  by  applying  the  refraction,  astronomic,  and  true  gravity 
corrections. 

Participation  in  the  redefinition  of  the  North  American 
Vertical  Datum  by  each  country  will  have  many  rewards.  Automation 
of  the  survey  observations  will  allow  the  network  to  be  upgraded 


by  eliminating  known  systematic  errors  in  the  net.  Possibly 
more  important,  is  the  future  use  of  the  data  by  scientists. 

As  new  leveling  data  are  observed,  vertical  crustal  motion  and 
rates  of  elevation  change  can  be  determined  by  comparison  with 
an  existing  data  bank  of  information  in  automated  format. 
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